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Abstract 

In a dense star, the Pauli exclusion principle functions as an enormous energy 
storage mechanism. Supersymmetry could provide a way to recapture this energy. If 
there is a transition to an exactly supersymmetric (susy) phase, the trapped energy 
can be released with consequences similar to gamma ray burst observations. Previous 
zeroth order calculations have been based on the behavior in a prototypical white 
dwarf of solar mass and earth radius (such as Sirius B) and have neglected density 
inhomogeneity. In this article we show that the effects of density inhomogeneity and of 
variations in masses and radii are substantial enough to encourage further exploration 
of the susy star model. In addition, the effects discussed here have possible applications 
to the growth of bubbles in other phase transition models in dense matter. 



1 Introduction 

In the past couple of years there has been a burst of theoretical activity discussing possible 
transitions between string vacua of differing amounts of supersymmetry [Q. These include 
a string theory study of a transition from a local minimum of positive vacuum energy to 
an exactly supersymmetric phase j2J. It is reasonable, therefore, to consider possible phe- 
nomenological signals for such transitions. We live in a broken-susy phase with, apparently, 
a positive vacuum energy density 
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leading to an acceleration in the expansion of the universe. The basic string theories, on 
the other hand, suggest a true vacuum of exact supersymmetry with massless ground state 
supermultiplets. In flat space these have zero vacuum energy. It is, therefore, interesting 
to consider whether we are living in a false vacuum which will ultimately decay to an exact 
susy ground state. False vacuum decays were treated in some generality many years ago jSj • 
In a false vacuum with energy density e it is expected that bubbles of true vacuum with 
radius r and surface tension S will have an effective potential given by a sum of a volume 
term and a surface term, 

Viyac) = -47r(r 3 e/3 - r 2 S) . (1.2) 

For small e the surface tension S can be treated as a constant. Bubbles will be constantly 
nucleated from the vacuum with a steeply falling distribution in initial radii. However, only 
those with initial radii greater than some critical radius 

R c (vac) = — (1.3) 
e 

will grow to effect a phase transition while smaller bubbles will be rapidly quenched. The 
probability per unit time per unit volume to produce a bubble of radius R c or greater and, 
therefore, to effect a phase transition to the true vacuum has been given [3| in the form 

r\ 2 P 

Ae -B(vac) (L4) 



dtd 3 r 

where, assuming a thin wall between the phases 
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In this picture, the fact that our broken susy world has existed so long is due to the smallness 
of e 3 relative to S 4 . In the physical vacuum, once a bubble of critical radius is nucleated 
the surface will rapidly accelerate and engulf the entire universe. The fact that this has not 
happened as yet suggests jlj that 

R c {vac) > R galaxy « 4.7 • 10 20 m (1.6) 

or 

S > 5.6 • 10 23 MeV/m 2 = 2 ■ 10~ 23 M o R E 2 . (1.7) 

In much of this article we use the solar mass, M = 1.2 ■ 10 60 MeV, and earth radius, 
Re = 6.38 • 10 6 m as convenient units. We also use natural units h — c—1. 

Reasonable expectations Ej exist that vacuum decay will be accelerated in dense 
matter. Heuristically, this can be seen by noting that, in the vacuum case, e is the energy 
advantage per unit volume of making a transition to the exact susy phase. One might expect, 
therefore, that the above equations will be modified in dense matter by replacing e by the 




Figure 1: The effective potential showing the false vacuum of broken susy and the true 
vacuum of exact susy. 



energy advantage per unit volume of trading the broken susy phase for the exact susy phase, 
i.e. 

e^e + Ap (1.8) 

where Ap is the ground state matter density in the broken susy phase minus the ground 
state matter density in the exact susy phase as shown in fig. ^ 

The difference Ap is the fermionic excitation energy density. The parameter controlling 
the exponential factor in the transition rate would then be 

27tt 2 S 4 , . 

b =wta^ ■ < L9 > 

The value of Ap in a white dwarf star is calculated as follows. In a degenerate electron 
gas of N electrons in a volume V the Fermi momentum is 



'3tt 2 N* 



1/3 



with, assuming equal numbers of neutrons and protons, 

Mn being the nucleon mass. The average kinetic energy is 

<E> -m = m(-l + 2 F 1 (-l/2,3/2;5/2;-p|/m 2 )) . (1.12) 
In the limit of zero electron mass which we will use for simplicity, this is 

<E>=^ . (1.13) 
4 

The kinetic energy density, which is equal to the difference in ground state energy densities 
between the broken susy state and the exact susy state, is then 

3p F N 1 /3vr 2 p\ 4/3 

We could double this estimate since a comparable contribution is expected from nuclear 
excitation energies but, for the present, we will neglect corrections by factors of order a few. 
From n~9l we have 

(L15) 

with 

1/4 

SM N . (1.16) 
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The phase transition probability per unit time per unit volume, Ae~ , increases rapidly 
with the density of the medium until the density becomes of order p at which point it 
saturates. For more dense media, the transition rate is proportional to the volume. 

The longevity of the universe (eq J1.7|) implies that 

p > 0.U0M o Re 3 (1.17) 
not far from the average density of the prototype white dwarf 

p WD = ^M Rf . (1.18) 

This suggests that the susy phase transition rate per unit volume could be appreciable for 
white dwarf stars but negligible for less dense objects and for denser objects of much smaller 
volume. The relative transition rates in a variety of astrophysical and terrestrial objects 
treated as of constant density are tabulated in table H One can see here that if p is close 
to its lower limit, the transition rate in white dwarfs is orders of magnitude greater than in 
the other considered bodies. In the future, it would be of interest to explore larger values 
of p for which the transition rates for white dwarf stars and neutron stars are comparable. 
However, since the density distributions in neutron stars and the radiative processes after 
the neutron to sneutron pair conversion are greatly different from the current calculations, 
these investigations are beyond the scope of the present paper. 
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Table 1: Masses, volumes, mean densities, and relative transition rates for a variety of 
physical bodies assuming p = 0.14. Masses and distances are given in units of solar mass 
and earth radius. 

In two recent articles |H] we have explored the possibility that such a transition in 
a dense star is the central engine of gamma ray bursts. The intense, collimated gamma 
radiation released in this way could provide the power to accelerate a macroscopic portion 
of the star to relativistic energies as in the cannonball model jH]- For reviews of the susy 
star idea, see [TUJUI]. In these articles we have presented zeroth order predictions based on 
a transition in a typical white dwarf star (e.g. Sirius B) neglecting density inhomogeneity 
effects. In the current paper, we proceed to incorporate these effects as well as consequences 
of existing variations in dense star mass and radii. 



In dense matter one would expect the critical radius to be 
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e + Ap(r) 
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In an inhomogeneous medium, a bubble of radius r will grow as long as 



r > R c (r) 



(1.20) 



Even ignoring density inhomogeneity within the star, the critical radius could be quite small 
inside the star but jump at the surface to a value that is much larger than the radius of the 
star. This has the effect of efficiently confining the susy phase to the interior of the star. 

In the exact susy phase the particle and sparticle have a common mass which we have 
assumed to be that of the particle in the broken phase. In order for the phase transition 
to proceed, it is a crucial assumption that the common mass is no greater than the particle 
mass. This assumption is, perhaps, supported by the fact that the ground state masses 
in the exactly supersymmetric string theories are zero. Also, susy models with radiative 
electroweak symmetry breaking (EWSB) relate the EWSB to the susy breakdown so that, 
in the absence of susy breaking, all particle masses vanish. There are no comparable physical 
models suggesting the opposite assumption, namely that the common mass in the exact susy 
phase is at higher energy. Nevertheless, this opposite assumption could also be considered 
but perhaps only briefly since, in this case, there is no exothermic phase transition unless, 
perhaps, the Fermi energy was greater than the mass difference. 

If a SUSY bubble forms in a dense material medium, the interior of the bubble will 
find itself greatly out of equilibrium since many of the fermionic constituents occupy high 
energy levels whereas scalar particles could all occupy ground state energy levels. In such 
a situation particle pairs will rapidly convert to sparticle pairs. Electrons, for example, will 
pair convert into scalar electrons (selectrons) thus evading the energy storing property of the 
Pauli exclusion principle: 



Note that, although the final state is a selectron pair and not a selectron-antiselectron pair, 
this process does not require R parity violation. 

Other types of phase transition models have also been proposed to explain gamma ray 
bursts. Among these are transitions to quark matter J2]> transitions to a color supercon- 
ducting state ^3], and transitions to mirror fermions ^3]. The calculations of the current 
paper on the behavior of a susy bubble can also be applied to hypothetical bubbles of these 
phases if they begin in the high density regions. 

The susy phase transition model applied to the typical white dwarf star neglecting density 
inhomogeneity predicts, correctly though roughly, in a relatively parameter-free way, the 
minimum duration of the burst (the light crossing time), the mean gamma ray energy (the 
mean electron kinetic energy), and the total burst energy (the total electron kinetic energy). 
In addition (and in distinction to the other phase transition models mentioned above), the 
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fact that the final state of the susy transition consists of scalar constituents predicts a 
significant amount of jet collimation due to Bose enhancement. 

In section II of this paper we review the density profile of white dwarf stars produced 
by the balance of inward gravitational pressure and outward electron degeneracy pressure. 
The electron momentum distribution differs greatly from that of a degenerate Fermi gas at 
uniform density. In section III we treat the classical free collapse time of an inhomogeneous 
susy star relieved of Pauli blocking, taking into account the full range of white dwarf stars. 
We also compute the variations in burst duration assuming that, in a dense star, the true 
vacuum bubble expands at the density dependent speed of sound which, in the high density 



2 Density and momentum space gradients in a white 
dwarf star 

Since Chandrasekhar it has been axiomatic in astrophysics that isolated stars below 
a mass of about 1.41 solar masses are stable against collapse due to electron exchange 
degeneracy. The density profile of such white dwarf stars is determined by hydrodynamic 
equilibrium between gravity and this outward degeneracy pressure initially augmented by 
thermal pressure. 

The degeneracy pressure of the electrons in a white dwarf takes a simple form in either 
the extreme relativistic or non-relativistic limits. In intermediate regimes, it can be written 
[To] in terms of the variable 



where \i e = A/Z = 2 and is the nucleon mass. The degeneracy pressure is proportional 
to the function 




Section IV is reserved for conclusions. 
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Specifically, 



P d = af(x) 
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with 
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The degeneracy pressure gradient is 
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Figure 2: The radius versus mass relation for a white dwarf at zero temperature determined 
by the balance of gravitational and electron degeneracy pressure. R and M are measured in 
units of earth radius and solar mass respectively. The prototype white dwarf, Sirius b, is 
indicated. 



In low temperature equilibrium this must balance the gravitational pressure gradient 



dP G p(r)G N M(r) 
dr r 2 

Here, M(r) is the mass within radius r and Gn is the gravitational constant. The resulting 
integro-differential equation for p can be solved by choosing an arbitrary starting value po 
at the center of the star and integrating outward until the density falls to zero, recording at 
each step the value of M(r). This defines the radius R of the star and the corresponding 
mass M(R) as a function of the peak (central) density. The resulting mass-radius relation 
[To] is shown in fig. El In standard astrophysics, all isolated stars with a mass of less than 
1.41M D will ultimately decrease in radius as they cool until they reach a point on the curve 
of fig. 121 at which point they become absolutely stable. In fig. |3] we plot the central and 
average density of white dwarfs at zero temperature as a function of the stellar mass. The 
central densities are more than an order of magnitude greater than the average densities. 
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Figure 3: The density at the stellar center (upper curve) and the average density (lower 
curve) in units of M R^ 3 plotted against the total mass of the star in units of solar mass 

The density distribution of these white dwarf stars as a function of distance from the 
center is a family of curves of which a representative seven are illustrated in fig. Previous 
work on the susy phase transition in dense stars 0IH] has ignored the strong density variation 
and relied on average densities only. The stellar density goes to zero at the surface of the star 



so, as can be seen from eq. 11.201 the susy bubble will stall at some distance from the surface 
creating a thin atmosphere of normal matter. The thickness of this shell is determined by the 
surface tension for which, at present, we know only the limit of eq ll.7l The above formulae 
imply, in this limit, a skin thickness of about 100/im for a star of solar mass and earth 
radius. The interior of the bubble constitutes a resonant cavity whose scalar constituents 
will continue to radiate gammas until the star radiates all its excitation energy or collapses 
under gravitational pressure. 




Figure 4: Densities as a function of distance from the center of seven representative dwarfs. 
The seven masses, in units of the solar mass, are 1.369, 1.330, 1.242, 1.072, 0.811, 0.512, and 
0.271. The dwarfs of higher mass have higher central densities and greater density gradients. 



In the inhomogeneous white dwarf, there is a local Fermi momentum given as in ll.101 
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The momentum distribution in the electron sea is no longer simply quadratic but takes the 
step function dependent form 
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The peak electron momentum and therefore the peak photon energy after the phase transition 
is given by eq. 12.71 evaluated at the stellar center, r = 0. These peak energies, tabulated in 



table 121 are much greater than that of white dwarfs treated as of constant density as in ref. 
[7] . Even higher peak energies, of course, can be found in stars closer to the Chandrasekhar 
limit since, there, the radius approaches zero and the density diverges. In this region there 
is a gradual approach to a neutron star structure. The gamma ray burst observations, for 
comparison, suggest a total burst energy of about 5 • 10 50 ergs and mean photon energies 
between 0.1 and 1.0 MeV. 

In table |2] we record various properties of these seven. 
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Pmax 
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1 


0.271 


1.92 


0.056 


0.009 


0.31 


0.013 


3.5- 


10 48 
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0.512 


1.44 


0.297 


0.041 


0.54 


0.034 


1.7- 


10 49 
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0.811 


1.0505 


1.572 


0.167 


0.93 


0.078 


6.1 ■ 


10 49 
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1.072 


0.7376 


8.321 


0.637 


1.63 


0.158 


1.7- 


10 50 
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1.242 


0.5002 


44.053 


2.366 


2.84 


0.291 


3.5- 


10 50 
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1.330 


0.3269 


233.238 


9.077 


4.94 


0.510 


6.6- 


10 50 
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1.369 


0.2062 


1234.875 


37.212 


8.63 


0.877 


1.2 • 


10 51 



Table 2: Masses, radii, central densities, average densities, maximum electron momentum 
, average electron kinetic energy, and total electron kinetic energy for seven representative 
white dwarfs at zero temperature. Units of mass and length are solar mass and earth radius 
respectively. The peak electron momentum is given in MeV/c, the average electron kinetic 
energy is given in MeV and the total electron kinetic energy is given in ergs. 



3 Factors affecting burst duration 

In the susy phase transition model there are several physical effects influencing the burst 
duration. These are 

1. The bubble growth time. This is the time it takes for a bubble nucleated in a high 
density region near the center to grow to the stellar surface. A lower limit to this time 
is the light crossing time but, more plausibly, in dense matter the bubble surface should 
expand at some rate comparable to the speed of sound in matter of that density. 

2. The light crossing time. After the bubble has engulfed the star, there could be an 
additional time required for light emitted on the far side of the star to cross the 
stellar diameter. At high density, these escaping photons might undergo a random 
walk leading to a time proportional to the square of the radius divided by a mean 
free path. However, it has been observed [IT] that the Landau-Pomeranchuk-Migdal 
effect [TS] would greatly increase the transparency of dense matter to gamma rays or, 
equivalently, the mean free path. 

3. The free collapse time. The conversion of fermions to bosons following a susy phase 
transition eliminates the degeneracy pressure and the star will undergo gravitational 



collapse. Classically, the star would collapse to a point in a matter of seconds but, once 
a radius a few times the Schwarzschild radius is achieved, general relativistic effects 
begin to dominate and will cause collapse to the Schwarzschild radius, as seen by a 
distant observer, to require an infinite amount of time during which photon emission 
will be progressively red-shifted leading to a certain amount of afterglow below the 
gamma ray spectrum. In ref. jH] it was found that the electrons higher in the Pauli 
tower have a higher probability of pair conversion implying an earlier pair conversion 
than that of lower energy electrons. This implies a natural progression to smaller 
frequencies as the burst progresses. The growth of the susy bubble into lower density 
regions has similar consequences. These effects by themselves have a much shorter 
time scale than observed afterglows and therefore are relevant only if the gravitational 
collapse and bubble cooling are greatly slowed by the following consideration. 

4. Radiation pressure. Radiation released by nucleons converting to scalar nucleons, 
which then drop into the nuclear ground state, is expected to rapidly expand the 
star and, thus, decrease the stellar density. This will increase the bubble growth time, 
the light crossing time, and the stellar collapse time. Density waves could be created 
in the wake of the initial blast and might lead to the rapid time variability observed 
in gamma ray bursts. Existing studies of susy bubble behavior, including the present 
article, do not include the effects of radiation pressure. 



We will study first the bubble growth time assuming this is governed approximately by the 
speed of sound in dense matter. The speed of sound at radius r depends on the local pressure 
and density and is given by 



v s (r) = ] jj- • (3.1) 



In the case of constant density this takes the simple form 



v s (r) = J^ 1 7rG N p(RZ-r*) . (3.2) 



Here 7 is the ratio of specific heats (5/3 for a monatomic gas), and R is the stellar radius 
at which the pressure vanishes, eq. 13.21 may be derived by considering the downward force 
exerted on a column of matter from radius r to the surface. The bubble growth time is then 

r R tiR 
r= / dr/v s (r) = -— . (3.3) 
Jo 2v s {(J) 



This time is about 2 s for our typical white dwarf but given the variations in density 
among the full sample of white dwarfs, the growth times based on the average densities as 
recorded in table El have a ratio of maximum to minimum of about 148. This is close to the 
observed ratio of short bursts but far from the observed ratio of about 10, 000 including long 
bursts and short bursts together. In any case we know from figure 0] that the densities are 
rapidly varying especially for the higher mass dwarfs. In the case of non-constant density the 



pressure at radius r remains governed by the differential equation of eq. 12.61 One begins with 
zero pressure at the stellar surface and integrates inward to find the pressure as a function 
of distance from the center. The local speed of sound is given by eq. 13. II subject to the limit 
cj in the high density limit. The bubble growth time is given by the left-most equality 
of eq. 13.31 . The stellar mass, radius, and bubble growth time are then given as a function 
of the central density. The probability per unit time of a susy phase transition in a star of 
radius R is 

— = 4nA r 2 dre ( "^' . (3.4) 

N dt Jo y ' 

Because of the exponential suppression at low density, p is a key parameter in identifying 
the gamma ray burst progenitor in the phase transition model. Although A is a free pa- 
rameter at this point, the dominant contributions to the bursts are likely to come from the 
largest bodies of density p or greater. The assumption here is that A is not strongly density 
dependent. 

It is interesting to explore the possibility that p is close to its lower limit from eq. 11.171 
This would imply, on the cosmological time scale, an imminent end to our current world 
of broken symmetry. If, on the other hand, p is very much greater than its lower limit, 
the bursts would come primarily from denser objects than most isolated white dwarfs, i.e. 
either neutron stars or larger bodies in the process of gravitational collapse. In this case, 
density enhancement through accretion such as in the collapsar model [T^], could play a role. 
However, since the Fermi momentum of the electron sea is proportional to p 1 / 3 as in eq. I1.1UI 
the peak photon energy potentially gives, in the susy star model, an upper limit on p. 

In the current paper we will explore the possibilities that p is related to its lower limit 
by a factor of 1,5, or 25. The rate of bursts as a function of total progenitor mass is given by 
multiplying 13 .41 bv the number of white dwarf stars of given radius or, equivalently, of given 
mass. We use the Sloan survey mass distribution given by [20J assuming that the relative 
numbers do not change greatly when extrapolated to zero temperature. The distribution is 
strongly peaked at 0.56 solar masses as shown in figEl 

The observed burst rate is about 

dN 711 

— « 5 • 10 - V" 9^ ■ ( 3 - 5 ) 

but the local rate is much lower suggesting that the fuel of gamma ray bursts has been 
largely exhausted and we are now primarily seeing bursts that happened long ago in more 
distant galaxies. The observed sample of bursts with identified redshifts seems to cluster 
around redshift 1. This interpretation is, perhaps, supported by the observed shortage of low 
luminosity dwarfs in our galaxy [21] and the, at first glance contradictory, apparent excess 
of dark objects of white dwarf mass in the galactic halo 



From the total energy release, table H3 might suggest a p value near or slightly above 
3.5M R~^ 3 . However, the investigation of beaming and other contributions to the energy 
release is at too early a stage to rule out larger values of p. Since the bubble growth time 
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Figure 5: Approximate distribution of observed white dwarf masses following ref. j2Hj 
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1 


.271 


1.92 


1.42 


12.3 


7.82 


3.5 • 10 48 
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.512 


1.44 


0.693 


5.80 


3.70 


1.7- 10 49 


18 


1.9-10-7 


4.4 • 10" 44 


3 


.811 


1.05 


0.348 


2.86 


1.82 


6.1 • 10 49 


0.78 


12 


7.9 • 10~ 8 
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1.07 


.737 


0.181 


1.47 


0.93 


1.7- 10 50 


0.46 


20 


31 


5 


1.24 


.500 


0.095 


0.76 


0.49 


3.5 ■ 10 50 


0.19 


11 


44 


6 


1.33 


.327 


0.049 


0.39 


0.25 


6.6 ■ 10 50 


0.062 


4.4 


24 


7 


1.37 


.206 


0.024 


0.19 


0.12 


1.2 ■ 10 51 


0.017 


1.4 


9.2 



Table 3: relative burst rates with given bubble growth times, constant density growth times, 
and classical collapse times, as well as total burst energies. The burst rates are given with 
separate arbitrary normalizations for three values of p. The peak of the probability distri- 
bution is underlined. 




is a monotonic function of dwarf mass, we can use the mass distribution of fig. |5] and the 
transition probability of eq. 13.41 to calculate the shape of the growth time distribution. This 
is shown in table. El for the three chosen values of p. It is encouraging that the growth times 
are in the range of the observed short burst durations but a critical test of the model must 
await the incorporation of the other factors influencing the burst duration as enumerated 
above. 

One of these factors is the collapse time of a star relieved of Pauli blocking. In [7] we 
have noted that the classical collapse time of a susy white dwarf is 



Although that paper considered only stars of constant density, the collapse time remains the 
same for inhomogeneous stars as long as one uses the average density. The classical collapse 
times of the representative sample of white dwarfs is given in table El However once the star 
approaches the Schwarzschild radius, general relativistic effects dilate the collapse time as 
seen by a distant observer. The approach to the Schwarzschild radius, r* = is given 

by 

r - r* = (r - r *)e- {t ~ to)c/r * . (3.7) 

During this time the star can still radiate although the photon energies will be redshifted 
from their emission energies according to the relation 

e -(t-t )c/2r* 

E = E em —j= ■ — . (3.8) 

J{l +e -(t-*o)c/r*) 

The time dependence of this component of the afterglow is independent of frequency. The 
time constant, c/r* for this contribution to the afterglow will be of order of ten microseconds 
for a star of near solar mass. With currently available techniques, it will be impossible 
to observe such a short afterglow. Another source of afterglow will be the emission from 
circumstellar material irradiated by the burst. This second source may be absent if the 
burst comes from the decay of an isolated star. Since, at present, afterglows have only been 
definitely observed for some of the long bursts, it might be interesting to consider a proposal 
where all of the short bursts and many of the long bursts originate in isolated stars while 
the bursts with extensive afterglows originate in stars with significant circumstellar material 
as in the collapsar model. A very massive star in the process of gravitational collapse will 
necessarily pass through stages of high fermion degeneracy where the possibility of a susy 
phase transition might become significant even though the time spent in these stages is not 
long. 



4 Conclusions 



We have considered the effects of mass variations and density inhomogeneities in white dwarf 
stars under the assumption that such stars experience a phase transition to the exact susy 
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ground state. It is clear that much work remains to be done. Nevertheless, it is encouraging 
that the basic assumption with few free parameters produces a gamma ray emission on the 
right time scale, with mean gamma energy and total energy release close to observations. The 
density inhomogeneities increase the expected peak photon energy by an order of magnitude 
over that of constant density stars of the same average density. 

The assumption that the surface tension of a susy bubble is independent of density needs 
to be examined and probably relaxed. Similarly, the strong density dependence of eqs. 11.411.91 
might be somewhat softened by corrections to the thin wall approximation of the vacuum 
decay studies [3]. Effects due to the radiation pressure at non-zero temperature need to 
be incorporated and may significantly affect the duration distribution. Radiation from the 
collapse of the Pauli tower in nuclei as well as the enhanced energy release from snuclear 
fusion and beta decay need to be studied. The resulting radiation pressure is expected to 
slow the final collapse of a star relieved of Pauli blocking. In addition, we hope to explore 
the possibility that the profound nuclear explosion at the stellar center due to the phase 
transition sets up density standing waves. This could be related to the spikey behavior of 
the observed gamma ray bursts. 
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